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An integrable eight state supersymmtric U model is proposed, which is a fermion model with 
correlated single-particle and pair hoppings as well as uncorrelated triple-particle hopping. It has an 
(?Z(3|1) supersymmetry and contains one symmetry-preserving free parameter. The model is solved 
and the Bethe ansatz equations are obtained. 
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Exactly solvable models of strongly correlated fermions have in recent years generated a great deal of attention, 
since they are believed to play a promising role in unraveling the mystery of high- Tc superconductivity (see, e.g. 
iQl). Several integrable correlated fermion systems have so far appeared in the literature. Most famous are the 
supersymmetric t-J model and the Hubbard model. Other integrable correlated electron systems of interest include 
the extended Hubbard model |0] , the supersymmetric U model proposed in |^ and extensively investigated in [^-^ , 
and its g-deformed version the model proposed in Q and its generalization jlOj . 

In this communication, we propose an eight state version of the supersymmetric U model. It is a supersymmetric 
fermion model with correlated single particle and pair hoppings as well as uncorrelated triple-particle hopping. We 
then solve the model by means of coordinate Bethe ansatz method and the Bethe ansatz equations are derived. 

Let (cj,a) denotes a fermionic creation (annihilation) operator which creates (annihilates) a fermion of species 

a — +, 0, — at site j. These operators satisfy the anti-commutation relations given by {cl^jCj^p} = SijSap, where 
i,j — l,2,---,L and a,(3 — +, 0, — . At a given lattice site j there are eight possible states: 

|0), c]^+|0), c]_o|0), c]_|0), 

4+4ol0>' 4+4-10)' 4o4-io)' 4+4o4-io)- (1) 

By rij^a = ^l^Cj^a we denote the number operator for the fermion of species a at site j. 

In the sequel we only consider periodic lattice of length L. The Hamiltonian for our new model reads 



^ 2{g\ 1) ^ (4"4/3Ci+i,/3Cj+i,a + h.c.) exp |-|(nj, ^ + nj+i,^)| 
- (g + l)(g + 2) (4+4o4-^.-+l,-^. + l^oC, + l,+ +h.C.) 
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'/ = -ln^, C-ln(.g + l)-iln.g(5 + 2), ^--In^. (3) 
5+1 2 5+2 

As will be seen, the supersymmetry algebra underlying this model is gl{3\l). Remarkably, the model still contains 
the parameter g as a free parameter without breaking the supersymmetry. Also this model is exactly solvable on the 
one dimensional periodic lattice, as is seen below. 

Our local hamiltonian Hij{g) does not act as graded permutation of the states (|^) at sites i and j. If one projects 
out any one species, then one singlet, two doublets and one triplet are projected out and the remaining four states 
are: one hole, two singlets and one doublet. In this case, the projected hamiltonian is nothing but that of the 
supersymmetrie U model Q with the parameter U (in the supersymmetric U model) related to the parameter g here 
via U = It is in this sense that the present model is an eight state version of the supersymmetric U model. 

To show that is gl{3\l) supersymmetric, we denote the generators of gl{3\l) by i?^, n,v = I, 2, 3, 4 with grading 
[1] = [2] = [3] = 0, [4] = 1. In a typical 8-dimensional representation of gl{3\l), the highest weight A = (0, 0, 0\g) itself 
of the representation depends on the free parameter g, thus giving rise to a one-parameter family of inequivalent irreps. 
Let {|a;)}|^i denote an orthonormal basis with |5), |6), |7) even (bosonic) and |2), |3), |4), |8) odd (fermionic). Then 
the simple generators and {-E^+i, £'/^~''^}Jl=i are 8x8 supermatrices of the form 
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The non-simple generators are obtained from the simple ones by using the defining (anti-)commutation relations of 
f/Z(3|l), which we omit. Further choose 

|1> = |0), |2)=4+|0), |3)=4o|0), |4)=4J0), 

|5)=4+4ol0>, |6)-4^.4J0), |7)=4ctjO), 18)^4+44 JO). (5) 

Then the verification that the hamiltonian H{g) commutes with all generators of gZ(3|l) is just a straightforward 
calculation. 

The model is exactly solvable by the Bethe ansatz, since the local hamiltonian Hi^i^i{g) is actually derived from a 
(?/(3|l)-invariant rational R-matrix [which satisfies the graded Yang-Baxter equation]. The R-matrix is given by [|llj 

p. ^ p u + 2g - iu + 2g){u + 2g + 2) ^ {u + 2g){u + 2g + 2){u + 2g + i) ^ 

' u-2g ' {u-2g){u-2g-2) ^ {u - 2g){u - 2g - 2){u - 2g - 4) ^' ^ ' 

where Pk, fc = 1, 2, 3, 4, are four projection operators: 

P,=j2\'^l){^ll F4-El*t>(*'.l, P2 = l{I + P)-Pu P, = \{I-P)-P,, (7) 

k=l k=l 

where P is the graded permutation operator and |^^.), 1^'^), k — 1, 2, • • • , 8 are given by 

|vl/}) = |l)®|l), |vI/i) = -L(|2)®|l) + |l)®|2)), 

1*1) = -L(|3) ® |1) + |1) ® |3)), 1*1) ^ -L(|4) ® |1) + |1) ® |4)), 
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1*1) = ^, [v/g+l(|5) ® |1) + |1) ® |5)) + V5(|2) ® |3) - |3) ® |2))], 

^2(25+1) 
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1*1) = ^, [Vg+l(|7) ® |1) + |1) ® |7)) + V5(|3) ® |4) - |4) ® |3))], 
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- ^^1_ [^{\2) ® |7) + |7) ® |2) + |5) ® |4) + |4) ® |5) ~ |3) ® |6) - |6) ® |3)) + ^F+^m <E> |1) + |1) ® |8))], 

= ^75=^[Vff + 2(|7) ® |2) - |2) ® |7) + |5) ® |4) - |4) ® |5) - |6) ® |3) + |3) ® |6)) + ® |8) - |8) ® |1))], 

*2) = ^W7^=^[V3+^(|8) ® |2) + |2) ® |8)) + ^^{\5) ® |6) - |6) ® |5))], 
^2(25 + 3) 
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^i) = ,,,, , ,J 7g+^(|8) ® |4) + |4) ® |8)) + ^iT2{\6) (g> |7) - |7) (g> |6))], 
^2(25 + 3) 

*5> - ^(-|8) ® |5) + |5) ® |8)), = -L(_|8) « |6) + |6) ® |8)), 

*7> = ^(-|8) ® |7) + |7) ® |8)), = |8) ® |8), (8) 

which are easily seen to be orthonormal, so that 

(*Il = (l*l»^ (*ll = (l*'.))^ fc=i,---,8, 

Here [|a;)] = for even (bosonic) \x) and 1 for odd (fermionic) \x). 

On the L-fold tensor product space we denote R{u)j j^i — ® R{u) (g) I^i^-o-'^) ^ and define the focal 

hamiltonian by 



= -4(2g + 1) (Pik,+i + (P4)„+i + 2.9P,,,+i. (10) 

u=o 5 + ^ 



Then by (p|), (ph and (|^), and after tedious but straightforward manipulation, one gets, up to a constant, 
Hjj+i{g) = 2(g+i) ^^j+i (ff)' [This identity also indicates that H{g) commutes with all generators of gl{3\l), since 
the rational R-matrix R{u) is a (7Z(3|1) invariant.] 

As mentioned above, the system is exactly solvable by means of Bethe ansatz technique. We assume the following 
wavef unction 



N 




(^Pqi > • • • , ^Pq„ )exp\i'^ kp^ Xj , (11) 



where Q is the permutation of the N particles such that 1 < xq-^ < ■ ■ ■ < xq^ < L. Denote Xq = {xq-^ < ■ ■ ■ < xq,^}. 
The coefficients ^qq^ . - ^qq^ (^Pqi j ' ' ' i ^Pq„ ) from regions other than Xq are connected with each other by elements 
of two-particle S-matrix: 

t'(fci) - 0[k2) + ic 

where operator V12 interchanges the species variables ai and 012 (ai, a2 = +, 0, — ), the rapidities 0{kj) are related to 
the single-particle quasi-momenta kj by 9{k) = i tan(-|) and the dependence on the system parameter g is incorporated 
in the parameter c = e'' — l = l/g. The periodicity condition for the system on the finite interval (0, L) results in the 
Bethe equations for the set of momenta kj : exp{ikjL) = Tj, j = 1, ■ ■ ■ , N, where 

'^j = 'S'ij+i(fci; %+i) • • • Sj^N{kj,kN)Sj^i{kj, ki) ■ ■ ■ Sjj-i{kj , kj^i), j = 1, - ■ ■ ,N. (13) 

The meaning of Tj is the scattering matrix of the j-th particle on the other {N — 1) particles. So now the problem is 
to diagonalize Tj to arrive at a system of scalar equations. It is easy to show that Tj = r(A = kj), where 

t(A) = tro [S'o,i(A - fci) • • • Som{X - k^)] (14) 
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is the transfer matrix of the inhomogeneous g^(3)-spin magnet of N sites. The commutativity of the transfer matrix 
for different values of the spectral parameter A implies that Tj, j = 1, ■ ■ ■ , N can be diagonalized simultaneously. The 
Bethe ansatz equations are written in terms of the rapidities 9j = 6{kj) and 
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The energy of the system in the state corresponding to the sets of solutions {6j} and {A^} is (up to an additive 
constant, which we drop) E = — 2^^-^^^ cos kj. 

To summarize, we have presented an integrable eight state version of the supersymmetric U model. It is a 5^(311) 
supersymmctric fcrmion model with generalized hoppings. We have solved the model by the coordinate Bethe ansatz 
method and derived the Bethe ansatz equations. There are many things remained to be done for this new model. 
One of them is to study physical properties, such as the phase diagram and the critical exponents, of the system. It is 
also interesting to incorporate integrable boundary conditions into the model and to investigate finite-size corrections 
of the boundary system. We hope to report results on those aspects in future publications. 
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